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results are compared with numerical simulations. 
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1. Introduction 

In the last two decades a large variety of noise induced phenomena have been discovered 
in far from equilibrium nonlinear systems. Among them there are stochastic resonance 
P, resonant activation [2j, noise enhanced stability and noise delayed decay [E]-[Ej- 

Recently it was predicted theoretically that the presence of additive noise may increase 
the average escape time from metastable states |H El CHI and from unstable states 
p3IH|. Moreover the noise can enhance the stability of an unstable fixed point in nonlinear 
discrete maps [3 E] ■ These resonance-like phenomena show a nonmonotonic behaviour 
of the average escape time as a function of the noise intensity. This means that by 
varying the noise intensity we can lengthen or shorten the lifetime of unstable states. In 
previous studies in fact was found that fluctuations can only accelerate the escape time 
from unstable states PSj-JSj- In the present work we study the phenomenon of Noise 
Delayed Decay (NDD) in more detail by analytical and numerical methods. Besides 
the Mean First Passage Time (MFPT), which quantifies only one of many possible 
characteristics of the decay process from unstable states, we consider the Nonlinear 
Relaxation Time (NLRT). The NLRT takes into account the inverse probability current 
directed from the stable state to the unstable one, which is neglected by the MFPT. 
The analytical calculations show that NDD effect for the NLRT is much greater than 
that for MFPT [Sj. In the present paper these analytical predictions are verified with 
numerical simulations. We show that the NLRT can be increased by the noise, when the 
noise intensity is varied in a very wide range. We consider the model of one-dimensional 
overdamped Brownian motion in the potential field of force 



Here x is the coordinate of the Brownian particle or the representative phase point 
denoting the state of the system, <&(x) is the potential describing the system itself, 
£(t) is the white Gaussian noise, with (£(£)) = 0, (£(£)£(£ + t)) = 2qS(r)/r], 2q/rj is 
the intensity of fluctuations, rj is the coefficient of equivalent viscosity, and q = kT 
is the energy temperature of fluctuations. The simplest case where the NDD appears 
is the parabolic potential. This potential profile was considered in theoretical papers 
concerning the decay from equilibrium unstable states fF2\ IT3].[T5]-[THj 



The initial state is an unstable equilibrium state, if Xq = 0, and it is an unstable 
nonequilibrium state if xq ^ 0. Various authors considered the decay in this system 
starting only from one point xq = 0. In this case, the deterministic decay time is 
infinity and the action of the noise decreases the decay time in accordance with known 
scaling laws ^ n ^ ne °t ner hand, the effect of the NDD always appears when 

we consider any xq ^ 0. In this case the decay times (both MFPT and NLRT) growth 
with noise, reach some maximum, and decrease. This effect was predicted theoretically 
in references [7J |H] • Before we consider the analytical expressions for MFPT and NLRT 




(1) 




(2) 
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let us describe in more detail their definitions. These definitions involve the decision 
interval restricted by one or two boundaries which specify the area of the unstable 
state. The MFPT is the average time that the Brownian particle stays within the 
decision interval before it crosses any of the absorbing boundaries for the first time. 
Due to the absorbing boundary conditions, the MFPT neglects that the particle may 
return into the decision interval after it has crossed the boundary once. On the other 
hand, the NLRT is an average time which takes into account this inverse flow directed 
inside the decision interval, when the particle comes back. The definition of MFPT is 
well known and it is a very useful time characteristic of decay processes, because the 
analytical expression for MFPT can be written in quadratures for an arbitrary potential 
profile (See e.g. references [IH]-|2II|)- Let P(t) be the probability that the particle is 
within the decision interval. Initially P(0) = 1. With time the particle escapes from 
the unstable state, therefore P(oo) = 0. If during the decay process the particle may 
cross the boundaries of the decision interval any number of times, the definition of the 
NLRT reads 



If the particle can cross the boundaries of the decision interval only once, the time © 
coincides with the MFPT. Otherwise the expression (j^J) takes into account the inverse 
probability current across the boundaries of the decision interval and differs from the 
MFPT. Therefore the MFPT is a particular case of the NLRT (|3*|) . namely, the case in 
which we neglect the inverse probability current. The analytical expressions for NLRT 
in quadratures for an arbitrary potential were obtained in [2T] . 

All numerical simulations shown in the paper were done using the algorithms 
described in [21], using the Heun algorithm for the stochastic integration, and the 
Ziggurath and a carry and subtract algorithm for the generation of the noise deviates. 
The integration time step used and the number of averages taken depend on the actual 
numerical experiment carried out, and on the set of parameters studied. Both have been 
changed throughout the numerical experiments aiming at optimizing the simulations: 
for example, when the noise intensity q was changed, the integration time step h was 
varied accordingly, to keep the quantity qh small, which is a prerequisite to have good 
convergence in the Heun integrator. 

2. Decay times for the symmetric potential 

For the symmetric potential profile $(— x) = the MFPT equals to (See e. g. 




(3) 



[231125!) 




(4) 



Let us consider the polynomial potential profile 

$(x) = -ax k /k, 



(5) 
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Figure 1. The MFPT normalized to the deterministic decay time versus the 
dimensionless temperature a = q/\&(L)\ for non-equilibrium states described by 
parabolic potential, for two initial states: xq = 0.6L and xq = 0.8L. Symbols are 
the results of numerical simulations and solid lines are the theoretical results obtained 
from equation (jSJ. 



where k = 2n is even. The exact expression for MFPT with potential (0) was 
obtained in ^2j) for a particular case, when the initial state is unstable equilibrium one 

Xq = 

where 2-^2(^1, a 2] &i> b 2 ; z) is a generalized hypergeometric function [2HII2ZI- It follows 
from (J3J, that for an arbitrary Xq (— L < x < L) the MFPT reads 

T(x , L; q) = T(0, L; q) - T(0, x ; q), (7) 

For example, if the potential is parabolic, then k = 2 and MFPT is 

2 



T(x ,L;q) = ^- 



x 2 02 F 2 1,1;-, 2, 



3 n ^^o) 1 



(8) 



2 ' 9 , 

When the noise is absent the MFPT coincides with the deterministic decay time In 
particular, for the parabolic potential the deterministic decay time is 

T d (x ,L) = (7] /a) ln(L/x ). 

The plots of the MFPTs normalized to for the parabolic potential and for different 
xq are shown in figure 1, where the symbols are the results of digital simulations and 
solid lines are the theoretical predictions. One can see that the NDD effect appears. The 
agreement between theory and numerical simulation is very good. To take into account 
the influence of the inverse probability current on the decay time we should consider 
NLRT. To obtain the NLRT we use the exact expression in quadratures obtained in 
references [HI for the symmetric potential 

r(x ,L;q) = T(x ,L;q) + Q(L;q), (9) 
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Figure 2. The normalized NLRT versus the dimensionless temperature a = q/\$>(L)\ 
for non-equilibrium states described by a parabolic potential, for two initial conditions, 
namely xq = 0.61/ and xo = 0.8L. Symbols are the results of numerical simulations 
and solid lines are the theoretical predictions. 



where T(x , q) is the MFPT © and 



Q(L;q) = 1 I" e^ v)/q dv ( 
a jl Jo 



~^ v)lq dv. 



For the polynomial potential (0) the NLRT correction ()10|) reads 



Q(L; q) 



-l/k 



-a 



i r 



1 



-a 



(10) 



q \k' 
where j(cx} a) is the incomplete gamma function j2jj and a = q/\<b(L)\. In particular, for 
the parabolic potential profile 0, when k = 2, the NLRT correction can be expressed 
in terms of error functions as 



e(L;q) = -j- (1 - Erf^; 



Erfi(a~ 



(12) 



The plots of NLRT are shown in figure 2, where symbols are the results of digital 
simulations and solid lines are the theoretical predictions. The NLRT increases with 
noise and displays the NDD effect. The enhancement of NLRT by noise is much greater 
than that of MFPT. This is because of the influence of the inverse probability current. 
Indeed, the NLRT correction (fTU|) is positive for any q > 0. 

Note that the NLRT may exceed its deterministic value for a great range of the 
noise temperature, e. g. for xq = 0.6L this range is < q < 6|$(L)| and for xq = 0.8L 
this range is more than one decade of dimensionless temperature. The value |$(£)| is 
the maximal height of the potential profile within the decision interval. Thus, the NDD 
effect takes place for the NLRT even for a strong noise intensity, while for the MFPT 
it appears only for weak noise when q < |$(L)|. When q > 3|$(L)| the asymptotical 
expression may be obtained for NLRT 



T(x ,L;q) 
k 



V 



l + m z )o-- L + 



aL k ' 2 

r(i) 



1 'i 



a 



k + l 



(13) 
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Figure 3. The normalized MFPT and NLRT versus temperature a = q/\$>(L)\ for 
non-equilibrium state described by the cubic potential = —x 3 /3, afor xq — 0.8L. 
Symbols are the results of numerical simulations and solid lines are the theoretical 
predictions. 



where m = xq/L. The asymptotic (jlHj) . normalized to the deterministic decay time 
T^(xo, L) is shown in figure 2 by dashed curve. It follows from equation ()13|) that when 
<r 3> 1, the NLRT does not depend on the initial conditions, i.e. on Xq. This is not the 
case for MFPT, which is 

T(x 0l L;q)^^- 2 (l-m")a-\ (14) 

for a ^> 1. We can not see this effect from figure 2 because the NLRT is normalized 
to different deterministic decay time. Therefore we may conclude that for strong noise 
intensity, the inverse probability current removes differences in decay times caused by 
initial conditions. 



3. Asymmetric potential 

Let us refer now to the asymmetric potential profile and consider the polynomial 
potential (jHJ) where k = 2n + 1 is odd. All the states xo > of this system are 
nonequilibrium. In reference ^2] the exact and approximate expressions for MFPT 
were investigated for xq < but not for xq > 0. While when xq > the effect of NDD 
appears. We consider the case x > 0. Let the decision interval be R : [— oo;L], then 
the MFPT reads 

T(x , L; q) = 1 [ L e* {v)lq f e-^ u)/q dudv. (15) 

and the NLRT, in accordance with [2*T] . is given by equation where 

e(L;q) = ^ e^ v)/q dv ■ f L e-^ u)/q du. (16) 
q Jl J-qo 

Using the expression for T(0, L; q) obtained in [T2] it is possible to write the analytic 
expression for MFPT when xq > 



T(x ,L;q) = — 
2q 



\ k k k q 
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Figure 4. Potential profile i|19fl . 
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(17) 



-$(L) -$(x Q ) 



q \ay \KJ \K q q ) 

The NLRT correction (|16|) for the polynomial potential (0) is 



2/Jfc 



-i 



(18) 



The plots of MFPT dTTJ) and NLRT Q for a cubic potential (fc = 3) and x = 0.8L are 
presented in figure 3, where the symbols are the results of numerical simulation. We 
see that the NDD takes place. The NLRT is always greater than the MFPT because 
the inverse probability current always increase the decay time. In general, despite the 
apparent difference between the expressions for decay time in antisymmetric potentials, 
the features of the times are similar. 



4. Potential with barrier 



Let us consider now a smooth potential profile with a barrier 

■? 

= -— + bx (19) 
3 

where a = 1, b = 1 and boundary is located at the point L = 5 (See the figure 4). 
Then, the local minimum of this profile is located at x — — 1, and the local maximum at 
x = 1. The barrier height is E = $(1) — $(—1) = 4/3. The expressions in quadratures 
for MFPT and NLRT for potential profile of this type were obtained in [HI They 
coincide with the corresponding expressions for the antisymmetrical potential (|9~|). (|15|). 
and (|16|). We will consider two cases of initial positions: Xq = 1.5 and Xq = 3. Both of 
these initial states are situated behind the potential barrier. The barrier height "seen" 
by the Brownian particle in its initial position xq is AE = $(1) — <&(xq) (see figure 4). 
Therefore AE/E < 1 for xo = 1.5 and AE/E > 1 for xq = 3. In references (HE] it was 
shown that the behaviours of MFPT and NLRT have strong dependence on the value 
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Figure 5. The normalized MFPT and NLRT versus the dimensionless temperature 
for potential profile with the barrier l|19|) and for non-equilibrium initial unstable state 
xq = 3 (AE > E). Symbols are the results of numerical simulations and solid lines are 
the theoretical predictions. 



AE/E. Indeed, if AE/E > 1, in the deterministic limit (i. e. when q — > 0) we obtain a 
finite value for the escape time 



By increasing the noise intensity the particles can go towards the well and be delayed 
there. As a consequence, the MFPT and the NLRT increase with q and go through 
a maximum for some q* > 0. This is the typical case of NDD phenomenon considered 
above (see figure 5). However, if AE/E < 1 (see figure 6) the situation becomes more 
complicated. It is evident that if there is no noise (£ = 0), the escape times (both 
MFPT and NLRT) are equal to the deterministic escape time (J2*U|) . At the same time, 
it follows from equation (fT3|) that in the limit of q —>■ 0, the MFPT and the NLRT go 
to infinity and for q = a singularity appears. This singularity was demonstrated for 
the first time in references jH |3] . The asymptotic expression for the MFPT as q — > 
and E > AE > can be obtained from equation (fTo|) : 



where x m i n is the coordinate of the local minimum. For both above cases the plots 
of MFPT and NLRT are shown in figures 5 and 6, where the symbols are the results 
of digital simulations and solid lines are the theoretical predictions. One can see, that 
the agreement between theory and the results of numerical simulations is good only for 
initial position x = 3 and there is strong difference when x — 1-5- In the last case, 
when q — > 0, the theoretical curve goes to infinity, while digital simulation gives the 
deterministic time. We can interpret this effect as follows: as explained in reference [3] 
the theoretical curve goes to infinity because even a very small noise intensity can 
eventually push the particle, initially located behind the barrier, back into the potential 
well. Then, the particle will be trapped there for a long time, because the well is very 




(20) 
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Figure 6. The normalized MFPT versus the dimensionless temperature a = q/E for 
potential profile with the barrier (|19f) and for non-equilibrium initial unstable state 
xq = 1.5 (AE < E). Symbols are the results of numerical simulations and solid lines 
are the theoretical predictions. 



deep. These trajectories of the Brownian particles lead to the singularity of the MFPT 
for q — > 0. On the other hand, the observation time in a digital simulation is finite and 
for some value of q it becomes smaller than the average escape time of a particle trapped 
in the well. Therefore the simulated escape time has the maximum at this point (See 
the region B in figure 6). Besides, the ensemble of particles in a numerical experiment 
is also finite. That is why for very small q, when the probability for a particle to be 
trapped in the well decreases exponentially to zero, we do not observe such particles in 
the simulation. As a consequence, the average decay time in simulation becomes equal 
to the deterministic one at q — > (See the region A in the figure 6). 



5. Conclusions 



We have derived the exact expressions for MFPT and NLRT in the case of 
polynomial potential profiles and for arbitrary unstable non-equilibrium initial positions. 
Expressions for MFPT obtained earlier (See e.g. [T2]) are valid only for the particular 
case when the initial position is the unstable equilibrium state and the action of noise 
always decreases the decay times. We have obtained nonmonotonic behaviour of the 
decay times as a function of the noise intensity and have shown that the inverse 
probability current can increase the NLRT with respect to the MFPT. Moreover the 
range in which we obtain the NDD effect is larger for the NRLT than for the MFPT. An 
important characteristic of the NDD effect is that we can both accelerate or slow down 
the decay of unstable nonequilibrium states by varying the intensity of fluctuations 
in a larger range for NRLT than for MFPT. The numerical simulations are in good 
agreement with the theory for most of potentials and initial conditions considered, 
except the specific case of singularity in the MFPT and the NLRT for q — ► 0. This 
singularity appears for the non-equilibrium states in a potential profile with a barrier 
when AE < E (figure 6) and it is also responsible for the enhancement of the stability of 
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fluctuating metastable states [HI [SI [XO] - In this case due to the limitation in the particles 
number and observation time in simulation, we cannot obtain the singularity (infinite 
decay time) numerically. Therefore one should expect that the same difference found 
between theoretical values of the average decay times and numerical simulations for 
AE < E and q — > will appear in experimental results when compared with theoretical 
ones. 
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